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ABSTRACT 
We study the quiet Sun magnetic fields using spectropolarimetric observations of the infrared and visible 
Fe I lines at 6301.5, 6302.5, 15648 and 15653 A. Magnetic field strengths and filling factors are inferred by 
the simultaneous fit of the observed Stokes profiles under the MISMA hypothesis. The observations cover 
an intra-network region at the solar disk center We analyze 2280 Stokes profiles whose polarization signals 
are above noise in the two spectral ranges, which correspond to 40% of the field of view. Most of these 
profiles can be reproduced only with a model atmosphere including 3 magnetic components with very different 
field strengths, which indicates the co-existence of kG and sub-kG fields in our 1"5 resolution elements. We 
measure an unsigned magnetic flux density of 9.6 G considering the full field of view. Half of the pixels present 
magnetic fields with mixed polarities in the resolution element. The fraction of mixed polarities increases as 
the polarization weakens. We compute the probability density function of finding each magnetic field strength. 
It has a significant contribution of kG field strengths, which concentrates most of the observed magnetic flux 
and energy. This kG contribution has a preferred magnetic polarity, while the polarity of the weak fields is 
balanced. 
Subject headings: Sun: magnetic fields - Sun: photosphere 



1. INTRODUCTION 

The quiet Sun can be defined as that part of the Sun far 
from magnetic activity. If we exclude plages and network, the 
remaining quiet area covers 90% of the surface of the Sun, 
and remains almost unchanged during the solar cycle (e.g., 
iHarvev-Angle 1993). Few decades ago, this internetwork 
(IN) was thought to be almost devoid of magnetic fields and 
consequen tly, its relevance for the solar magnetism was not 
obvious ( Livingston & Harvev 1975; Smithson"1975V How- 
ever, with the present instrumentation, when the spatial reso- 
lution is aro und 1", the quiet Sun turns out t o be full o f mag- 
netic fields (^Grossmann-Doerth et al."1 996t iLin & Rimmeld 
[1399; Sanchez Almeida & Lites 2000), with an unsigned flux 
density measured with techniques based on the Zeeman effect 
of the order of 10 G (e.g., Sanchez Almeida & Lites 2000l 
JLites 2002). The analysis of observations with a spatial 
resolution twice bett er yiel ds an unsigned flux twice larger 
(iDominguez Cerden a et alJ l20 03a)'. There are reasons to 
believe that these resolutions are not enough to completely 
resolve the magnetic structures (see Sanchez Almeida et aU 
l2003b. and references therein). Recent measurements based 
on the Hanle e ffect find an unsigned flux of, at least, 60 G 
llTruiillo Bueno etaDl2004t iBommier et al. 2005), which is 
much higher than the magnetic flux observed with Zeeman 
techniques. Considering the fraction of solar surface corre- 
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sponding to the quiet Sun, the flux measured in the IN is com- 
parable to (or larger than) that coming from all active regions 
even in the maximum of activity. Thus, the quiet Sun may 
play a significant role in the physical processes governing 
the solar global magnetism (e.g ., Unno 1959; Stenflo 198^ 
'Sanchez Almeida'2003a';' Schriiver & Title 2003). 

The topology and the scale of variation of the quiet Sun 
magnetic fields are not yet clear There is evidence for the 
existence of mixed polarities, at least when the spa tial resolu- 
tion is I" (Sanchez Almeida et al. 1996; Sigwarth et al.1 19991 
I Sanchez Alm eida & Lites 2000; Lites 2002; Khomenko et alj 
12003b . The small-scale magnetic structures in the quiet Sun 
seem to be smaller than 50-1 00 km (see iKhomenko et alJ 
l2003t iDominguez Cerdena et al...2003b.) . This size assumes 
one single structure in the resolution element, however, it 
is to be expected that the dynamical influence of the gran- 
ulation on the magnetic fields produces a distribution of 
many smaller structures. This fragmentati on shows up in 
MHD simulations of turbulent dynamos (iCattan ed 119991 
Emonet & Cattaneo 2001), and 3 mag neto-convection (e.g., 
Istein & Nordlund 2002; V6gie?'200?, 'Vog ler et alj i2005i) . 
Recently, Cameron & Galloway (2005) studied various mech- 
anisms proposed to produce intense magnetic concentrations 
and limited the radii of the magnetic flux tubes to less than 
10 k m, which is appro ximately the magnetic diffusion length 
scale ('Schussler"1986). This scenario was already proposed 
by Sanchez Almeida et al. ( 1996) with an atmosphere formed 
by optically thin magnetic structures with diverse proper- 
ties (see also .Sanchez Almeida. .1998.) . Such a scheme is 
caUed MISMA (Micro-Structured Magnetized Atmosphere). 
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An inversion procedure based on it has been able to fit the 
whole variety of asymmetri c Stokes profiles^ produced by the 
quiet Sun magnetic fields ("Sanchez Al meida & LitesI 120001: 
[Spcas-Navarro & Sanchez Almeida 200^^! 

Observational studies based on infrared (IR) spectral 
lines find magn etic fie ld strengt hs in the range below 
1 kG (sub-kG; iLin & R immele [T999t rXhomen ko et alJ 
120031) . while works using lines in the visible spectral 
range t end to find kG field s iGro ssmann-Doerth et aD 
ISigwarth et al." "1999"; Sanchez Almeida & Litep 
Socas-Navarro & Sanchez Almeida 2002 

2003a b; Lites & Socas-Navarra 



TABLE 1 

Atomic parameters of the lines under study. 



1996 
200C 

Dommguez Cerdena et al. 
2004i) . Kilogauss fields are also inferred from the presence 
of a cons iderable number of G-band bright points in th e 
quiet Sun ("Sanchez Almeida et al."2004';'de Wiin et al."2005|). 
[Sdnchez Almeida & Lites (2000) put forward a solution for 
this apparent contradiction: the co-existence in the resolution 
element of a distribution of field strengths going from zero 
to kG. They show how to reproduce the observed bias with 
the vertical gradient of the field strength induced by the 
hydrostatic equilibrium in the atmosphere. Another possi- 
bihty with the same consequences is a h orizontal gradient 
dSocas-Navarro & Sanchez Almeidall2003l) . The presence of 
a range of field strengths from sub-kG to kG fields bias the 
measurements based on IR Stokes profiles to sub-kG fields. 
Such a particular behavior is due to the high sensitivity of 
the IR lines to magnetic fields. The Fe i A 15 648 A (used by 
ILin & Rimmelell999tlKhomenko et alJ2003[) is magnetically 
saturated^ at 400 G while the visible pair Fe I A6301.5 A, 
Fe I A6302.5 A (e.g., used by Sanchez Almeida & Lites 2000; 
iDominguez Cerdena et ali .2003^) saturates at 1.5 kG. The 
IR polarization signals coming from the strongest fields are 
spread in a wide range of wavelengths and they produce little 
signal. On the contrary, weak fields saturate the IR line and a 
significant Stokes V signal is built up in a narrow wavelength 
range. Weak field and strong field components appear in the 
same wavelengths in the visible lines so that their relative 
contribution is proportional to their magnetic flux. The weak 
field component can dominate over the strong component in 
the IR Stokes V profiles, while the strong fields dominate over 
the weak ones in the visible profiles. Thus, an incomplete 
analysis of the IR profiles leads to weak fields, whereas the 
same analysis of the visible lines shows strong fields. In 
order to test this conjecture, Sanchez Almeida et al. (2003a) 
presented simultaneous spectropolarimetric observations in 
the visible and IR ranges. An independent Milne-Eddington 
inversion of the visible and IR lines suggests the co-existence 
of sub-kG and kG fields. However, a more sophisticated and 
simultaneous inversion is desirable to provide further support 
to the conjecture of Sanchez Almeida & LitesI ( l2000l) and 
ISocas-Navarro & Slnchez Alm eida ( 2003^ 

This paper presents such simultaneous inversions under the 
MISMA approximation of the visible and the infrared Stokes 
profiles used in Sanchez Almeida et al. (2003a). As we dis- 
cuss above, it seems to be appropriate for reproducing the 

^ We use the standard Stokes parameters to define the polarization state: 
/ for the intensity, Q and U for the two orthonormal states of lineal polariza- 
tion and V for the circular polarization. The Stokes profiles of a spectral line 
give the variation with wavelength of these four parameters. If the spectral 
line is formed in a homogeneous atmosphere, the Stokes profiles are expected 
to be symmetric or antisymmetric (e.g., Landi DeslTnnocenti 1992). 

' The Stokes V signals no longer have a linear relationship with the mag- 
netic field strength. The level of polarization stays constant while the splitting 
of the Stokes V lobes grows in proportion to the field strength. 
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complex Stokes profiles produced by the quiet Sun. Section|2l 
summarizes the observations and the data reduction. Some 
properties of interest of the magnetograms obtained in the two 
spectral ranges are shown in § |3l The strategies for the data 
analysis and the inversion are explained in §|4| We present the 
results in §|5lwhose consistency is analyzed in §|6l A final §0 
discusses the main conclusions. 

2. OBSERVATIONS AND DATA REDUCTION 



We re-analyze the spectra used bv 'Sanch ez Almeida et alJ 
(2003a) whose properties are only sketched in the original 
reference. For the sake of comprehensiveness, the data and 
reduction are detailed here. The part of the reduction in which 
we bring the visible and IR spectra to equal spatial resolution 
is new. 

We use two telescopes for the simultaneous visible and IR 
observation"*. The observation was performed at the Span- 
ish Observatorio del Teide (Tenerife, Spain). The IR spec- 
tra were gathere d with the Tenerife Infrared Polarimeter (TIP, 
ICoUadosifT999t iMartinez PiUet et alJfT999l) operated at the 
German Vacuum Tower Telescope (VTT). Spectrograms of 
the full Stokes vector (/, Q, U, V) were obtained in the IR iron 
lines Fe l A 15648 A (Lande factor gL=3) and Fe I A 15652 A 
(gL=1.53). The visible observations were carried out with the 
French-Italian telescope THEMIS (Telescope Heliographique 
pour r Etude du Magnetisme et des Instabilites Solaires) using 
the spectropolar imetric mode MTR (Mein & Rayrole 1983 
iRavrole & Meini 119931) . The lines analyzed are the pair 
Fe I A6301.5 A (gL=1.66) and Fe i A6302.5 A (gL=2.5). The 
atomic parameters of the four spectral lines are listed in Ta- 
ble ^ In order to obtain spectra in a two-dimensional field of 
view (FOV), the solar surface was scanned from East to West 
with a step size of 0."5. The slit width was set to 0"5 in both 
telescopes and the scan consisted in 60 positions. To improve 
the signal to noise ratio, the integration time was set to 30 s 
per position. 

This work analyzes data obtained on August 10 2002 in a 
quiet Sun region close to disk center (// = 0.94). We choose 
a region without magnetic activity as judged with the help of 
Ca II K and Ha slit-jaw images from video cameras at the 
VTT. A video link between both telescopes providing the slit- 
jaw images from VTT at THEMIS was helpful to point both 

'' Nowadays, the observation is possible using a sin gle telescope with 
TIP+POLIS (VTT, Tenerife, see lKhomenko et alJl200.% and wi th SPINOR 
(DST, Sacramento Peak Observatorv. see lSocas-NavaiTO et all20 0Sl. 
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Fig. 1 . — Intensity images in the continuum of the IR lines before (left) and after (middle) smearing, and continuum image of the visible lines (right). The 
rectangles point out structures that can be easily identified in both the IR and the visible images. Axes are in arcsec from the lower left comer. 



telescopes to the same target. Since the absolute pointing of 
the telescopes was not accurate enough, and the slit-jaw im- 
ages of the granulation did not allow us to recognize the re- 
gion pointed by the two telescopes, we used another differen- 
tial method for accurate pointing. First, both telescopes were 
pointed to a characteristic small structure which could eas- 
ily be identified, e.g., a small pore. The pointing of the VTT 
was moved to a quiet Sun region, and then THEMIS pointing 
was displaced by the same distance. As we explain below, the 
scans provided by the VTT and THEMIS are post-processed 
for fine-tuning co-alignment. This final step removes residual 
pointing errors plus some additional effects (e.g., differential 
refraction between the visible and the IR). The scan under 
study here has an offset between telescopes of only 1" (see 
below). Since we scanned with a step size of 0"5 and with an 
exposure time of 30 s per step, the time lag between IR and 
visible data of the same point corresponds to only 1 minute. 

The data reduction of the two data sets was carried out in 
a similar way. A flat-field correction was applied to the spec- 
tra. In order to get the Stokes parameters, we demodulate the 
measurements with the calibration matrix, and finally we cor- 
rect for seeing-induced crosstalk by combining the two beams 
of each polarimeter. The crosstalk produced by instrumental 
polarization (IP) of the telescopes themselves was considered 
too. Since THEMIS is free of IP, the visible data do not need 
correction. Due to the configuration of mirrors in the VTT, 
the IR data are affected by an important crosstalk between the 
different Stokes parameters. It was calibrated with a model 
Mueller matrix set using a linear po larizer mounted at the en- 
trance window of the telescope (see ICollado^ 2003) . 

Since the purpose of this work is the simultaneous inver- 
sion of visible and IR lines, one needs a consistent scale of 
wavelengths in both spectral ranges. We carried out an abso- 
lute wavelength calibration so that the lines are unshifted in 
this scale when the solar plasma is at rest. The central wave- 
length of Fe I A6301.5 A and Fe I A6302.5 A was considered 
to be the zero-crossing of the Stokes V profiles correspond- 
ing to the signals in the network points of our FOV Such 
a reference gives a n acc uracy of 200 ms"' (Solanki 1986)). 
ISolanki & Stenflol (119861) find that the zero-crossing wave- 
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length depends on the spectral resolution. Using their sim- 
ulations, we check that our resolution of 21.9 mA per pixel 
does not threaten the absolute accuracy given above. The 
strongest polarization signals were chosen to trace the net- 
work (V > 0.021c with Ic the intensity of the continuum near 
the lines). Weaker signals tend to be red-shifted. Since this 
method has not been tested for the IR lines, their absolute 
wavelengths need another way of calibration. We computed 
the mean line profile from the data and then corrected for the 
expected convective blue-shift. We average the Stokes / pro- 
file over the FOV. The minima of the lines were fitted with 
fourth order polynomials whose minima were taken as the 
central wavelengths of the spectral lines. These wavelengths 
corrected for the convective blue-shifts provide the reference 
for an atmosphere at rest. Table |2l presents the convective 
blue-shift for Fe I A 15648 A and Fe I A 15652 A obtained from 
a two-com ponent invers ion of many line profiles in the quiet 
Sun ( Borre ro et al.l2003l) and fro m line synthesi s from numer- 
ical simulations of granulation jAsplund et alj|200 0. private 
communication). We cannot asses the accuracy of our abso- 
lute wavelength calibration based on them, however, the fact 
that the two numerical simulations provide the same convec- 
tive blue-shift suggests the consistency of the approach. 
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Fig. 2. — Magnetograms from Fe I A15648 A (left) and Fe I A6302.5 A (right). White and black distinguish the two polarities. Both magnetograms are scaled 
to flux densities |Beffl < 20 G. The distance between small tickmarks is 1". 



One of the most important steps in the data reduction is the 
spatial co-alignment of the visible and IR scans. The aim of 
this process is to ensure that the scans are not only taken from 
the same region on the Sun but also that each individual pixel 
of the visible and IR scans come from the same point of the IN 
region. The correct alignment must be guaranteed for a mean- 
ingful analysis of the data, i.e., for the simultaneous interpre- 
tation of both data sets. In order to align the scans we displace 
and re-scale them allowing for five degrees of freedom. First, 
we correct for differences in the spatial scales along the X and 
Y directions. We define X as the direction of scanning and Y 
as the direction along the slit. The step width of the scans 
was slightly different in both telescopes (less than 1 %), and 
the spatial scales along the slit were computed from the focal 
length of the telescopes. Their (relative) agreement with re- 
spect to each other were checked using reference images of a 
sunspot. Then, we found the shift in the X and Y directions 
by cross-correlation of the visible and IR continuum images. 
It corresponds to the offset of 1" mentioned above. Finally 
both slits were not exactly oriented with the same angle with 
respect to the Sun. There was a small rotation of 1° which 
we also determine by cross-correlation. Once the correspon- 
dence between the visible and IR spatial coordinates was set, 
we interpolate the visible spectra to obtain the spectra at the 
positions of the pixels in the IR scan. 

It is important to have the same spatial resolution in both 
data sets. However, the angular resolution of the observa- 
tions at the VTT was better than at THEMIS due to the use 
of a correlation tracker (Ballesteros et al. 1996). To solve this 
problem, the IR data were degraded to the resolution of the 
visible data by convolution with a Gaussian. The contrast 
of the IR continuum image depends on the angular resolu- 
tion. We set the width of the Gaussian, and so the angular 



resolution of the IR map, by comparison with the synthetic 
images from numerical simulatio ns of magneto-convection 
JVogleril2003t iVoder et alj l200?). Specifically, we tune the 
IR contrast so that the ratio between the visible and the IR 
contrasts agrees with the ratio obtained from numerical simu- 
lations (0.48). Such theoretical ratio was found to be almost 
independent of the angular resol ution and in excel l ent ag ree- 
ment with the value measured bv'Khom enko et aP J2005h us- 
ing a single telescope to obtain simultaneous visible and IR 
data. The central panel in Figure[2shows the continuum im- 
age in the IR after smearing. The contrast of the IR image is 
0.9%. The right panel of Figure^contains the continuum in- 
tensity image of the visible data (contrast 1.9%). The squares 
point out regions easy to identify in both images showing the 
goodness of the alignment. 

After all data reduction, including smearing of the IR 
data and Fourier filtering of the visible data, the level of 
noise in polarization for Stokes V is (2-3)xl0~'*/c for the 
IR data and 5xlO~''/c for the visible data. A more 

detailed description of the telescopes and t h e data reduc- 
tion can be found in iDominguez Cerdefia (12004 1). De- 
scriptions of the instruments can be foun d in ColladoJ 
(1999) and Schlichenmaier & Collados' ('2002') for TIP, and in 
Lopez Ariste et al. (2000) and Bommier &J^a.yra|fl, ( 2002 ) for 
THEMIS. 

3. MAGNETOGRAMS 

Figure |2] shows the magnetograms calculated from the 
Stokes V profiles of Fe I A6302.5 A and Fe I A15648 A. 
The longitudinal magnetic flux density or B^s was computed 
with the magnetograph equation und er the weak field approx- 
imation using the same method as Domin guez Cerdena et alJ 
(2003a) . The magnetograms thus obtained are equivalent to 
the routine magnetograms (e.g., that obtained with the MDI 
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instrument on board of SOHO^), if we exclude the velocity 
fields. Figure 121 displays both magnetograms with the same 
scale of magnetic flux densities. At first sight, both images 
are similar with the exception that the signals found in visible 
are larger than the signals found in IR. Indeed, this is the first 
important difference. Since only the IN fields are of interest 
here, the network region pointed by a circle in the lower left 
part of our FOV was discarded from the analysis by rejecting 
pixels with flux densities larger than 20 G in the visible and 
10 G in the IR. The noise levels turned out to be 2 G in the 
visible and 1 G in the IR. The method to estimate the noise is 
detailed in Dominguez Cerdeiia et al. (2003b). The unsigned 
flux density obtained from pixels in the FOV with flux density 
below the upper limit and above noise was. 



B, 



effi 



3.0G,forFelA6302.5 A, 
1.5G,forFelA15648A. 



(1) 



In both cases the pixels with signals cover 45% of the total 
FOV (the rest of the pixels were set to zero in the calculation 
of |Beff !)■ Under the same conditions of noise and spatial res- 
olution, the measurements in the visible detect more flux than 
in the IR. This suggests that visible and IR lines trace different 
magnetic structures. The difference of flux is conspicuous, 
and one may be tempted to ascribe it to the approximations 
used for Beff. However, more elabor ated techniques show 
the same be havior (ME inversions in .Sanchez Almeida et al.l 
l2003atlDorn inguez Cerdena 200l and §5). 

Some magnetic flux concentrations are very similar in both 
magnetograms (squares in the upper part of Figure |2j but 
some others differ (see the other two squares). There are 
structures having one polarity in the IR magnetogram and the 
opposite in the visible one. This effect is evident even in the 
network area with large polarization signals (see the circle in 
Fig. 13. Such a result can be considered as an additional evi- 
dence of the different sensitivity to magnetic fields of the vis- 
ible and the IR lines. 

The magnetic fields in the visible magnetogram (Fig. |2] 
right) exhibit a preferred polarity while the IR magnetogram 
(Fig. 121 left) looks well balanced. The signed flux was mea- 
sured to be. 



BeS = 



-2.0G, forFelA6302.5A, 
-O.IG, forFelA15648 A. 



(2) 



This imbalanced flux is similar to that mea- 
sured in other visible observations (e.g., Lites 2002; 
iDominguez Cerdena et al.i .2003^) . Yet, .Khomenko et al.i 
J2003h found almost no signed flux using the same IR lines 
(they observed 0.2 G). 

All these results indicate that visible and IR 

lines may trace dif ferent magnetic structures^ as 

conjectured bv iSan chez Al meida & Lites (i2000l) : 
ISocas-Navarro & Sanchez AlmeidaU2003l). 

4. DATA ANALYSIS 

In order to infer the magnetic field strengths from the ob- 
servations, we carried out an inversion of the Stokes pro- 
files under the MISMA hypothesis. Micro-Structured Mag- 
netized Atmospheres ( MISMAs) \y ere introduced for the 
first time by Sanchez Al meida et alJ ([1996), with the objec- 
tive of reproducing the asymmetric Stokes profiles observed 
in most solar structures. Unresolved optically thin micro- 
structures naturally render asymmetric Stokes profiles. The 

' See|http://sohowww.nascom.nasa.gov/| 



MISMA scenario turns out to be an appropriate scheme 
of analysis since the observed quiet Sun Stokes profiles 
are very asymmetric, and it has been shown to work with 
the asymmetries of Fe I A630 1.5 A and Fe I A6302.5 A 
(ISanchez Almeida & Lites 2000). In this section we give a 
description of the procedure followed to invert the data. 

4.1. PCA classification 

We only analyze those Stokes V profiles clearly above 
noise. A threshold of three times the noise level was cho- 
sen in order to separate clear polarization signals from noisy 
profiles. Pixels with profiles whose peak polarization is above 
the threshold in both the IR and the visible data cover 40% of 
the FOV. Thi s data set incl udes som e 22 8 pixe ls. 

Foflowing ISanchez Alm eida & Litesl dJOOOl § 3.2) flie 
Stokes V profiles in the selected data set were classified us- 
ing a principal component analysis algorithm (PCA). The aim 
is twofold. On the one hand, it reduces the number of rep- 
resentative profiles and, on the other hand, it provides low 
noise profiles to perform inversions used as initialization for 
the inversion of individual Stokes profiles. The classification 
is carried out simultaneously for the visible and the IR Stokes 

V profiles. The profiles of the four lines were scaled to the 
blue lobe of the Fe I A6302.5 A Stokes V profile. The PCA 
method gives 38 classes. All profiles with large polarization 
signals were included in one single class to avoid contami- 
nating the other classes with network profiles. The Stokes 

V threshold was chosen to be 1%/^ in the visible, and this 
special class includes all profiles from the network patch in 
the bottom left of our map (see Fig. |2j. The large number 
of classes compared with those obtained for visible lines by 
Sanchez Almeida & Lites (2000, 10 classes) or for IR lines 
by Kh omenko et al.l (2003, 8 classes) is due to the combi- 
nation of the different asymmetries found in the Stokes V 
profiles of the two pairs of lines. Note that, in general, the 
classes found in this work do not correspond to those obtained 
by Sanchez Almeida & Lites (2000) and by Khomenko et aU 
(2003). 

We represent each class with the average of all profiles 
corresponding to the class. The individual profiles are not 
scaled to carry out the average. We only multiply them by 
the sign of the blue lobe of Fe I A6302.5 A. Some of these 
classes are shown in Figures |3] and @] (the dots). There is 
a clear improvement in the signal to noise ratio (S/N, de- 
fined as the maximum of the unsigned Stokes V profile di- 
vided by the noise) from the individual profiles to those rep- 
resenting the class. The individual profiles have a typical S/N 
between 3 and 10, whereas the two first classes have a S/N 
of 100. The IR profiles from classes and 4 have tails ex- 
tending far out from their extremes, as those predicted by 
ISocas-Nava rro & Sanchez Almeida' ("2003). These tails come 
from the contribution of kG field strengths. Such profiles are 
similar to class in lKhomenko et alJ (I2003L Fig. 1). The first 
6 classes include 42% of the total number of profiles. One 
conspicuous fact is that some classes (e.g., classes 1 1 and 20) 
show one polarity in the visible and the opposite in the IR. 
Such a result was expected from the magnetograms in Fig- 
ure|2l Approximately 20% of the profiles show the effect (see 
also Sanchez Almeida et al. 2003 
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4.2. MISMA inversion 

The MISMA inversion code is described by 
ISanchez Almeidal ( Il997h . We will follow a strategy similar 
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Fig. 3. — Some representative classes of observed Stokes V profiles showing the large variety of asymmetries produced by the quiet Sun magnetic fields. The 
circles represent the observed profiles and the solid lines the fits with MISMA model atmospheres. The header on top of each row displays the class number, the 
percentage of profiles belonging to the class, and the magnetic field of the three MISMA components at the base of the photosphere, sorted in order of decreasing 
mass. The polarity of the magnetic fields is indicated by plus and minus signs. 



to lSanc hez Alm eida & Lit es ('200Cf). Here we summarize the 
main characteristics of the underlying model atmospheres: 

1 . There are four different components, one non-magnetic 
and three magnetic. Not e that the inversion of 
visible lines performed by Sanchez Almeida & LitesI 
(l2000h require two magnetic components. As 



ISanchez Almeida et alJ J19961 §4.2) discuss, it is 
impossible to reproduce even the mildest observed 
asymmetries with only one component. Here we 
need one additional component to characterize the 
weak magnetic field streng t h revealed by the IR 
lines (iLin i 
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Fig. 4. — Same as Figure|3]for 6 additional classes. 



iSanchez Almeida et al J 120033*). The components are 
embedded in a background producing unpolarized stray 
light. 

2. The magnetic fields are vertical, along the line of 
sight. Since the observed Stokes Q and U signals 
are rarely above the noise level, it is reasonable to 
assume the inclination of the magnetic field to be 



0° or 180°. The inversions were performed only 
for Stokes / and Stokes V. This approach allows 
us to compute longitudinal magnetic fields, however, 
the observations are compatible with inclined fields 
producing a Unear polarization hidden by noise (see 
ISanchez A lmeida & Lites"2000'. § 4.7). The impHca- 
tions of this assumption are discussed in S 16.21 



Dommguez Cerdena, Sanchez Almeida & Rneer 




-0,2 0.0 0.2 -0.2 0.0 0.2 

A-6301.499 [A] X-6302.492 [A] 



-0.5 0.0 0.5 
X-1564B.515 [A] 



-0.5 0.0 0.5 
A-15652.874 [A] 



Fig. 5. — Example of a fit of the Stokes / profiles of the four lines. The circles represent the observed profiles, and the solid lines are the fits. Note the two 
strong telluric blends in the red wings of Fe I A6302.5 A and Fe I A15648 A. 
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Fig. 6. — MISMA model from the inversion of class 0. The six pan- 
els show the variation with height of different physical parameters of inter- 
est. From left to right and top to bottom: the continuum optical depth at 
5000 A, the temperature, the magnetic field strength, the occupation frac- 
tion, the density, and the velocity along the field lines. The parameters are 
displayed for the four components except the optical depth: the sohd hne 
corresponds to the non-magnetic component, and the dotted, the dashed, and 
the dot-dashed lines to the three magnetic components. Note that the spec- 
tral lines are formed between 100 km and 300 km, so that the properties of 
the atmosphere outside this range are not constrained by the observations and 
therefore are not reliable. 



3. The temperature of the components is forced to be the 
same. Since the magnetic structures are optically thin 
the radiative exchange between components is very ef- 
ficient, washing out temperature differences in seconds. 

4. The four components are in lateral pressure balance. 



5. Motions are only allowed along magnetic field lines. 

The temperatures are retrieved at four fixed heights in the 
atmosphere. The parameters characterizing each compo- 
nent are set at the bottom of the model atmosphere. They 
are the magnetic field strength, the occupation fraction, the 
bulk velocity, and the broadening parameter. The stratifi- 
cation of these parameters is obtained under the assumption 
of hydrostatic equilibrium and horizontal pressure balance. 
Three parameters account for the stray light, and another one 
parametrizes both the solar macroturbulence and the instru- 
mental broadening due to limited spectral resolution. De- 
spite the fact that the spectral resolution in the visible and 
the IR differ by 30 %, one can use a single macroturbulence 
for the two wavelength domains because the macroturbulence 
inferred from the inversions is completely dominated by true 
solar motions rather than instrumental smearing. The to- 
tal number of free parameters amounts to 25. They are used 
to fit some 400 observables (each wavelength position of the 
Stokes / and V profiles of the four lines). 

An example of the fit to the intensity profiles is shown 
in Figure |5] The observed profiles (circles) exhibit line 
blends. These are a telluric line close to the red wing of 
Fe I A6302.5 A, a solar Fe I line at 15647.3 A, and an uniden- 
tified possibly telluric line in the red wing of Fe 1 A 15648 A. 
To avoid the contamination by these blends, the Stokes / pro- 
files at the wavelengths of the blends are excluded from the 
fit. There are also two OH bands of solar origin close to 
Fe I A 15652 A, however, these bands are only visible in cool 
solar regions such as sunspots. The intensity profiles in Fig- 
ure|5]are those from class 0. 

All models have three magnetic components. In order to 
distinguish between them, we denote them as first, second 
and third components. The first component is that carrying 
the largest mass, as was defined by Sanchez Almeida & Litej 
(2000), and it is followed by the second and the third. If 
the three components would have the same magnetic field 
strength, the first would be the one producing the largest po- 
larization signal. 

Figure |6l displays the model MISMA of class (see the 
Stokes V in Fig. O. The different physical parameters of 
the non-magnetic, the first, the second, and the third mag- 
netic components are plotted with solid, dotted, dashed and 
dot-dashed lines, respectively. The first magnetic component 
presents a weak magnetic field strength (^100 G at the base 
of the photosphere), while the other two have strong field 
strengths (larger than 1 kG). The density of the first mag- 
netic component is the same as for the non-magnetic com- 
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ponent (Fig.|6lbottom-left panel) whereas the other two com- 
ponents have much lower density. The reduction is forced by 
the coupling between field strength and density produced by 
the lateral pressure equilibrium. The higher the magnetic field 
the lower the gas pressure and, for the same temperature, the 
lower the density. 

5. RESULTS 

Figures |3] and 13 show the fits produced by the MISMA in- 
version of 12 representative classes. All kinds of observed 
profiles are properly reproduced with three magnetic compo- 
nents per resolution element. 

From now on, all results correspond to a single height in 
the atmosphere. We select the base of the quiet Sun photo- 
sphere, defined as the height where the total pressure equals to 
1.3 X 10^ dyncm"^. Invoking lateral pressure balance among 
the individual model MISMAs, this definition sets a single 
geometrical height for all the inversions, and it corresponds 
to the height where the continuum optic al depth is one in th e 
ID quiet Sun model atmospheres (e.g.. 'Maltbv et al.lll986h . 
The study is focused on two particularly important physical 
parameters, namely, the magnetic field strength and the oc- 
cupation fraction. The latter corresponds to a volume filling 
factor 

The magnetic field strength inferred from the MISMA in- 
version of each class is also given in Figures |3] and 0] The 
magnetic field of the components are denoted as Bi, B2 and 
B3, in order of decreasing mass. Very often the inversions 
require one component with a sub-kG field strength and two 
with strong kG fields. Such a behavior is typical of the most 
abundant classes (e.g., classes 0, 1, 4). The Fe I A15648 A 
Stokes y profiles of these clas ses are similar to that expected 
from the synthesis by Socas-Navarro & Sanchez Akneidal 
(|^03). The lobes of the profiles have a splitting equivalent to 
that produced by a weak field, however, they exhibit extended 
wings characteristic of strong fields. Some other classes have 
all components with fields in the range clearly below 1 kG 
(e.g., classes 5, 19, 20). In this case, there are no extended 
wings in the profiles of Fe I A15648 A. Finally, note the pres- 
ence of very asymmetric profiles. The IR profiles tend to have 
a higher degree of asymmetry. Some 28% of the analyzed 
pixels present IR profiles with three lobes ( e.g., class 5). We 
found a fraction similar to that observed bv lKhomenko et al] 
(I2OO3. 30%). There is also a small fraction of very asymmet- 
ric IR profiles with four lobes (3%, classes 14 and 30), while 
there are only a few visible profiles with three lobes (3%, e.g., 
class 30). 

The model atmospheres obtained from the classes were 
used to initialize the inversion of the individual Stokes V pro- 
files belonging to the class. From this moment on we refer 
only to results obtained from the inversion of the 2280 indi- 
vidual profiles. 

5.1. Magnetic field strength and magnetic fiux 

The histograms of the distribution of magnetic field strength 
obtained from the inversion of the individual profiles are 
displayed in Figure They show the percentage of pix- 
els within bins of 50 G for the different components. The 
solid, the dashed, and the dotted lines distinguish the first, 
the second, and the third magnetic components, respec- 
tively. The results from the network patch were omitted 
(see § |3}- The first component tends to have much weaker 
fields than the other tw o and it has a shape sim ilar to the 
distribution obtained by iKhomenko et al.l J2003h using only 




Fig. 7. — Histograms of the magnetic field strength! at the base of the pho- 
tosphere, determined from inversion of the individual Stokes / and V profiles. 
The first, the second and the third components are plotted with solid, dotted 
and dashed lines, respectively. 



IR lines. The mean field obtained for this component is 
290 G. The second and the third components show the same 
behavior as the histograms obtained with visible lines by 
Sanchez Almeida & Lites (2000). They have mean fields of 
970 and 1 100 G respectively. Despite the similarity between 
the magnetic field strengths of the second and the third com- 
ponents, they are well separated by the inversion code. There 
are large differences of Doppler shift between them; 85 % of 
the models are separated by more than 500 m s"' at the base of 
the photosphere, and this difference increases substanti ally i n 
the heights of formation of the spectral lines. (See also 8^ 16.31 ) 

The top panels of Figure |8l show the maps of field strength 
for the three magnetic components. The magnetic field of 
the first, the second and the third components are displayed 
from left to right. The network point (bottom left) is almost 
the only feature with strong kG fields in the first component. 
Most of the pixels show at least one component with strong 
fields. More specifically, 92% of the analyzed pixels possess 
at least one field larger than 500 G while the fraction with 
field strengths above 1 kG is 74%. 

The map with the filling factor of the magnetic fields (the 
sum of occupation fractions of the three components) is dis- 
played in the bottom left panel of Figure|8l It goes from 10% 
in network points to 0.5%, and the mean value is approxi- 
mately 4%. If we consider the whole FOV (all values in the 
pixels not analyzed are set to zero), and we exclude the net- 
work profiles, the measured magnetic fields fill only 1.5% of 
the IN region. 

The bottom right panel of Figure |8] gives the map of the 
total magnetic flux density computed as the product of the 
occupation fraction times the magnetic field strength. Note 
that the flux density is given in G (or Mxcm"^), thus it can 
be directly compared with other measurements with differ- 
ent pixel size or FOV. Most parts of the analyzed surface 
present a flux density larger than 25 G (sum of the three com- 
ponents). Considering the whole map, we measure a total 
flux densit y of 9.6 G which is similar to the value obtained 
bv Socas-Navarro & Sanchez Almeida (2002) or Lites ( 2002) 
with better spatial resolution. 
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Fig . 8 . — The three upper panels (a,b,c) show maps of the magnetic field strength at the base of the photosphere of the first, the second and the third components, 
respectively. The bottom left panel (d) gives the occupation fraction (filling factor) of the three components all together in a logarithmic scale. The total unsigned 
magnetic flux is shown in the bottom right panel (e). It is saturated at 30 G. Black regions correspond to pixels with polarization signals too low to be inverted. 
The separation between minor tickmarks is one arcsec. 
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Fig. 9. — Mean magnetic field strength of all components as a function of 
the magnetic flux density. The shaded area gives the standard deviation above 
and below the mean value. 



Figure |9] shows the variation of the magnetic field strength 
of all components with the magnetic flux density. The shaded 
area represents the standard deviation above and below the 
mean value. A difference between the two standard devia- 
tions reveals the asymmetry between the distribution of field 
strength above and below the mean value. There is a strong 
tendency for the field strengths to weaken when the magnetic 
flux decreases. Sanchez Almeida & Lites (2000) obtained a 
less pronounced dependence, but their lowest detectable flux 
density was 10 G while in this work it decreases down to 1 G. 

5.2. Mixed polarities 

50% of the analyzed pixels can be fitted only with mixed 
polarities in the model MISMA. Since we fit 40% of the total 
pixels, this fraction corresponds to 20% of the FOV. This large 
amount of pixels includes all profiles with three or four lobes 
in the IR and/or visible (30%), and also those profiles with a 
clear polarity in the visible Stokes profiles and the opposite in 
the IR profiles (20%). We find that the fields of those pixels 
are weaker than the rest. Figure^] shows the histograms of 
the field strength as in Figure 0but only for models including 
mixed polarities. All components show fields weaker than in 
Figure \J\ From the first, the second and the third magnetic 
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Fig. 10. — Same as FigureQbut only for those pixels with mixed polarities 
in the resolution element. 
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Fig. 1 1 . — Percentage of pixels with mixed polarities versus the Stokes V 
signal of Fe I A6302.5 A. 



components we retrieve a mean field of 250, 700 and 880 G, 
respectively. 

Figure^2gives the percentage of mixed polarities as a func- 
tion of the Stokes V signal of Fe 1 A6302.5 A. The weaker the 
polarization signal the larger the fraction of mixed polarities. 
This behavi or was already noted by^S anchez Almeida & Lites 
( l200Qi) and iSocas-Navar ro & Sanchez Akneidai (i200 2) but 
they obtained a maximum fraction of 25%. The fraction in 
Figureiniflattens to 65% for Stokes V lower than 0.002. 

5.3. Distribution of field strengths 

We compute the probability density function (PDF) of the 
magnetic field strength from the magnetic field and occupa- 
tion fractions procured by the MISMA inversions of the in- 
dividual profiles. The PDF gives the probability of finding 
a specific magnetic field strength in the observed IN region. 
The range of magnetic fields from to 2000 G was divided 
into bins of 50 G. Each individual intrinsic field B is assigned 
to one bin. The probability to find the magnetic field of a 
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Fig. 12. — Probability density function of the magnetic field strength at the 
base of the photosphere as inferred from Zeeman signals. 



single bin is the sum of the occupation fractions of all pixels 
with B belonging to this bin. We normalize the PDF so that 
its integral provides the filling factor of magnetic structures in 
the IN region assuming that the pixels without signals have no 
magnetic field strength. 

The PDF of the magnetic field strength condenses informa- 
tion on the quiet Sun magnetism that can be used for compar- 
isons with other observations and numerical simulations. The 
first moment (B) is connected with the magnetic flux density 
and the second moment provides the magnetic energy density 
(B2)/87r: 



(B): 



Jo 



B PDF(B) dB, 



B^ PDF(B) dB. 



(3) 



(4) 



Figure E] shows the PDF of our IN region. Once again, 
we have avoided those pixels belonging to the small network 
patch in the FOV. The shape of the PDF is to be expected from 
the histogram of Figure Most of the filling factor is con- 
centrated in the weak fields, while most of the magnetic flux 
comes from the strong fields. Table |3l summarizes the prop- 
erties of the PDF. It shows the filling factor, the flux density, 
the fraction of flux, and the fraction of energy in four differ- 
ent ranges of magnetic field strengths, i.e., the full range, the 
weak fields (B < 500 G), the strong fields (B > 500 G), and the 
kG fields (B > 1000 G). Three quarters of the total observed 
flux and 90% of the magnetic energy is in the kG fields, de- 
spite they represent only 30% of the total filling factor 

There is an accumulation of field strength between 1600 
and 1800 G. It may be due to the combined effect of a mag- 
netic intensification mechanism plus the existence of an upper 
limit that the magnetic field strength cannot exceed. This limit 
of 1800 G corresponds to a magnetic pressure equal to the gas 
pressure of the non-magnetic plasma at the base of the photo- 
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TABLE 3 
Filling factor, magnetic flux density and magnetic energy density of the PDF 

IN FigureFT^ 



Range 


filling factor 


flux 


density 


energy density 


fraction of flux 


fraction of energy 




[%] 




[G] 


[ergcm"^^] 


[%] 


[%] 


total 


1.50 




9.6 


500 


100 


100 


B<500 G 


0.87 




1.2 


10 


10 


2 


B>5()0 G 


0.63 




8.4 


490 


90 


98 


B>1000 G 


0.42 




7.0 


450 


75 


90 



sphere. The magnetic intensification mechanism would try to 
pile up magnetic fields right before the upper limit. 

5.4. Signed flux 

When we average over the FOV, the signed magnetic flux 
does not cancel completely. The flux density, taking into 
account the direction of the magnetic field, turns out to be 
-3.5 G. This net flux wa s expected from the results in § |3l 
(equation 13). iLitesI (2002) already observed the presence of 
an imbalanced flux in almost all the observation he studied. 
He measured an imbalance (defined as the signed flux divided 
by the unsigned flux) between -0.20 and +0.48, and found a 
direct correlation between the imbalance of the IN fields and 
the imbalance of the network fields around them. Our im- 
balance is -0.36 and our network region is negative almost 
everywhere. Figure^jshows the PDF of the signed magnetic 
field strength, i.e., the magnetic field strength times the sign 
of the vertical component of the magnetic field. There is a 
large contribution of strong fields in the negative flux part of 
the PDF, while the positive part has a continuous fall off from 
weak to strong fields. The weak fields are well balanced while 
the strong fields have a preference toward negative polarity, 
i.e., the same polarity as the network patch. For fields weaker 
than 500 G the imbalance is only -0.09, and for field strengths 
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Fig. 13. — PDF of the magnetic field strength times the sign of the polarity 
of the flux. The vertical dotted line points out B = G. 



Stronger than 1 kG the imbalance is -0.5."Socas-Na varro et alJ 
(2004) already observed the preference of the strong fields to- 
wards one polarity, and Figure[0]shows a clear observational 
evidence of such an effect. 

6. CRITICAL ASSESSMENT OF SOME ASSUMPTIONS AND 
RESULTS 

Some of the properties of the observational material used in 
the work, and the assumptions of the inversion, are not ideal 
but dictated by feasibility. In order to strengthen the validity 
of the results, we devote this section to analyze the assump- 
tions and results that may be more critical. 

6.1. Time lag between IR and visible data sets 

Our visible and IR spectra average the IN properties in spa- 
tial scales of l."5 arcsec and temporal scales of 30 sec (§|3. 
The time lag between visible and IR spectra would be of con- 
cern only if this spatio-temporal average changes significantly 
during 1 min. Fortunately, large changes of the average prop- 
erties are not expected for a number of reasons. First, the in- 
tegration time of the individual measurements smears out the 
short time-scale variations of the IN fields, and only variations 
in the range between 30 sec and 1 min are of concern. Sec- 
ond, low spatial resolution observations of IN magnetic fields 
show how man y structures rem ain stabl e for tens of minutes, 
and longer (e.g.JZhang e t al.ll99&.Dominguez Cerdena et alJ 
l200 3b). Obviously, our time lag does not represent a prob- 
lem for these structures. In view of the results coming from 
high resolution observations and numerical simulations of 
magneto-convection, this time stability cannot not be inter- 
preted as the stability of individual magnetic structures but 
the stability of the average properties. For example, a strong 
granular downdraft lives longer than the granules and, dur- 
ing its lifetime, it continuously advects magnetized plasma 
towards a specific point of the solar surface (Rast 2003). Ob- 
served with low resolution, one would detect a magnetic patch 
that remains in place for a long period of time, despite the 
fact that it is formed by many small-scale short-lived mag- 
netic structures that are continuously advected and engulfed 
by the downdraft. Finally, for the average magnetic properties 
to present significant variations during the time lag, the typi- 
cal velocities of the magnetized plasma would have to be un- 
physically high. Consider the horizontal velocities required to 
sweep out during the time lag all magnetic elements existing 
in a resolution element. In this case the magnetic structures in 
the same resolution element of the visible and IR maps will be 
completely different, questioning any interpretation based on 
their simultaneity. The required velocities would have to be 
of the order of the speed needed to cross a resolution element 
during the time lag, i.e., 18 km s~'. Such speed is too high 
to be common. It is supersonic, and much larger than the ob- 
served vertical velocities. Obviously, the lack of simultaneity 
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becomes more problematic as the spatio-temporal resolution 
of the observation improves. 

From our point of view, the spatial and temporal integration 
of the measurements represents a problem more serious than 
the time lag and the small spatial resolution differences that 
may remain between the IR and the visible data sets,. The 
polarization detected in each resolution element results from 
adding up contributions from many magnetic structures hav- 
ing very different properties. Our inversions provide only a 
crude representation of what may be happening in the real 
Sun. This problem, however, is not specific of our work but 
affects all previous observational works on the IN magnetic 
fields. In the case of the MISMA inversions the problem is 
ameliorated because a sketchy model MISMA actually repre- 
sents many different underlying atmospheres, including very 
complex ones. All those atmospheres having the same aver- 
age properties produce the same spectrum, and these average 
properties are the ones that the inversions retrieve (see the dis- 
cussion in Sanchez Almeida & Lites 2000, § 5). 

6.2. The inclination of the magnetic fields 

The inversions were carried out under the assumption of 
longitudinal magnetic fields, which minimizes the number of 
free parameters and simplifies the inversions. This assump- 
tion, however, biases the inferred PDF(B). Fortunately, the 
importance of the bias can be estimated and it turns out to be 
small. If a magnetic field has an inclination 7 with respect 
to our vertical line-of-sight, the Stokes V spectrum remains 
approximately as in the case of vertical field but scaled with 
the factor cos 7. This is a general property of the radiative 
transfer equation, and is n ot restrict ed to weak fields on plane- 
parallel atmospheres (see Sanchez A lmeida & Truiillo Buenol 
fl999. § 3.1). Consequently, if an observed Stokes V profile is 
reproduced with the wrong inclination, all the parameters of 
the inversion would remain the same except for the occupa- 
tion fraction 01(E). The occupation fraction provides a global 
scaling factor for the Stokes V profiles, and so it absorbs the 
full bias. Our PDF(B) is biased because it is based on occu- 
pation fractions underestimated by a factor | cos 7(B) |, 



n{B) 



PDF(B) =A?-' ^ |cos7,(B)|a,(B), 



(5) 



/■=! 



with the subscript / denoting one of the n{B) pixels with a field 
strength B, and the symbol A^ standing for a normalization 
constant. Obviously, PDF(B) sets a lower limit to the true 
PDF*(B), 

n{B) 

PDF* (B) = A^"' ^ aiiB) > PDF(B) , (6) 

(=1 
since |cos7,(B)| < 1. On the other hand, it is possible 
to set an upper limit by estimating the maximum inclina- 
tion allowed by the noise in Stokes Q and U. Following, 
ISanchez Almeida & Lites ( 2000, § 4.7), we synthesize the lin- 
ear polarization produced by our model atmospheres when the 
magnetic field of the three components is inclined. For each 
model MISMA we determine the largest inclination 7,, pro- 
ducing undetectable linear polarization, i.e., below the obser- 
vational noise. If we use this inclination to correct the occupa- 
tion fractions, one can define a new PDF„(B) which represents 
an upper limit to the true PDF*{B), 

n(B) . . 

PDF„(B) = A^-' V I °' ,J| «/(g) > PDF*(B), (7) 



|cos7„/(B)| 



where the inequality follows from the fact that 
|cos7,|/|cos7a,| > 1 for all /. We find that the unsigned flux 
corresponding to PDF„(B) is only 25% larger than the flux of 
PDF(B), indicating that the assumption of vertical fields have 
a very limited effect on the PDF estimated in the paper. 

6.3. Need for three magnetic components 

Due to the similarity between the magnetic field strengths 
of the second and the third magnetic components, one may 
wonder whether they are really needed to reproduce the 
line shapes. Moreover, if they are not needed, the exis- 
tence of mixed polarities may be an artifact since the sig- 
nals of the two components cancel our producing no ob- 
servable residual. As we stress in § 14.21 the three compo- 
nents are needed. Here we elaborate on the reasons. Re- 
producing the Stokes V asymmetries of the visible lines re- 
quires two magnetic components, even for the mildest asym- 
metries characteristic of the network. One single magnetic 
component produces Stokes V area and the peak asymme- 
tries of the same magnitude, in contr adiction with the observa- 
tions ( Sanchez Almeida et,^ljl22a p. 545 ). Somehow unex- 
pectedly, Sanchez Almeid a & LitesI OOOOh find that these two 
components are enough to explain not only the mildest asym- 
metries, but the whole range of observations including ex- 
treme profiles with an even number of lobes. These two com- 
ponents have similar kG magnetic field strengths, but very 
different velocities (Sanchez Almeida & Lites 2000, Figs. 12 
and 14). The large difference of velocities is demanded by the 
observed Stokes V asymmetries, and it is the reason why the 
inversion code is able to separate th em. Here we analyze pro- 
files similar to those reproduced by Sanchez Almeida & LitesI 
(I^POO) and, consequently, a minimum of two magnetic com- 
ponents is needed. In addition, we also want to reproduce 
the IR lines which, according to the literature, indicate the 
existence of weak fields in the quiet Sun. A third weak mag- 
netic field component was needed too. This argument does 
not imply that all the observed profiles necessarily require 
three components. It implies that the typical cases need three 
components (regular Fe I A6301.5 A andFe I A6302.5 A, with 
Fe I A15648 A and Fe I A15652 A showing weak fields). It 
also implies that that typical Stokes V profiles have enough 
information to constrain the three components. Finally, it im- 
plies that if a single MISMA scenario is going to be used to 
reproduce all the asymmetries, it must contain three compo- 
nents. 

6.4. The inferred magnetic field strengths 

One of the main conclusions of this work is the co-existence 
of kG and sub-kG fields in resolution elements of l."5 as 
inferred from the simultaneous inversion of IR and visible 
Stokes profiles. 

The use of visible lines to measure kG magnetic field 
strength has been criticized. The analysis of extremely noisy 
Stokes V profiles of Fe I A6301.5 A and Fe I A6302.5 A 
may induce a bias towards kG (see Bellot Rubio & Colladqj 
I2OO3I) . and the Stokes V profiles from the visible lines are 
the major tracers of the kG magnetic fields in our inversion. 
However, the kG that we infer are not produced by such a 
bias for a number of reasons. The classes have much less 
noise than the individual profiles, and yet the inversion yields 
as much kG of the individual profiles. The classes that yield 
at least one kG magnetic component represent 75% of the 
total pixels. They include the two main classes with a S/N 
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Fig. 14. — PDF of the magnetic field strength computed from the inversion 
of the IR profiles alone. The solid line is the PDF inferred from the inversion 
under the MISMA hypothesis, and the dashed line the PDF from the Milne- 
Eddington inversion. 



larger than 100. Another 10 classes have S/N larger than 
30, which ensures inversions wi thout the bias pointed out by 
IBellot Rubio & CoUadosI (12003! ). 

Another way to test the consistency of our results is to in- 
vert only the IR lines to see whether it is possible to infer 
kG fields in the quiet Sun without visible lines . Th is inver- 
sion was performed in the same way as in § 14.21 but only 
for Fe 1 A 15648 A and Fe 1 A 15652 A with MISMA mod- 
els including only 2 magnetic components (see more details 
in Dominguez Cerdeiia 2004). The result of such analysis was 
a major magnetic component similar to the first component in 
Figure [7] and a second with weak and strong fields. A field 
stronger than 1 kG was found in 50% of the inverted pix- 
els (those with the strongest IR signals). The strong fields 
were retrieved from the extended wing of the Stokes V pro- 
files. Figure^^gives the PDF obtained from these inversions. 
The strong field contribution is not as large as in Figure E] 
but it still carries most of the magnetic energy and flux. The 
presence of fewer strong fields as compared to the full inver- 
sions including visible lines is due to the incompleteness of 
the information given by the IR lines, since most of the con- 
tributions from strong fields remains hidden by noise. The 
signed flux ha s a b ehavior similar to that of the full inversion 
described in § 15.41 It is more important for the fields stronger 
than 1 kG while the weak fields are balanced. Figure [T4l also 
shows the PDF obtained from the Milne-Eddington inversion 
of the IR lines (the dashed line; see Sanchez Almeida et a^ 
l2003ah . This inversion was performed with a single magnetic 
component and shows only the weak field component. 

7. DISCUSSION AND CONCLUSIONS 

A new observational technique was used to study the quiet 
Sun magnetism. Simultaneous visible (Fe I A6301.5 A 
and Fe 1 A6302.5 A) and infrared (Fe 1 A15648 A and 
Fe 1 A 15652 A) polarimetric observations were analyzed un- 



der the assumption of Micro-Structured Magnetized Atmo- 
sphere (MISMA). We generate model MlSMAs able to re- 
produce all kinds of asymmetries presented in the observed 
Stokes y profiles, even when there are clear differences in the 
asymmetries of the visible and the IR Stokes profiles. The 
model atmospheres have three magnetic components. 

The simultaneous inversion of the visible and IR pairs re- 
veals the co-existence of kG and sub-kG fields in the low 
photosphere. We have observed this combination of weak 
and strong fields in 74% of the pixels und er st udy, which 
means 30% of the pixels of the FOV (see § 15. 1> . This was 
already observed by Sanchez Almeida et al. (2003a) from the 
separate inversion of the IR and visible lines, however, this 
work presents the first evidence of the co-existence of kG 
and sub-kG fields from a simultaneous inversion of lines in 
the two spectral ranges. From the three magnetic compo- 
nents, the one having most of the mass includes most of 
the weak fields. The two other components tend to show 
strong kG fields. The preference of the first component to 
have the weakest fields is due to the coupling between field 
strength and plasma density. The weaker the field the larger 
the density so that the component of weakest field tends to 
have the largest mass. The histogram of field strengths of 
these minor components (Fig. |7) looks sim ilar to that ob- 
tained by Sanchez Almeida & Lites| ( 120001) . while the first 
component has a distribution similar to th at obtained with the 
analysis of IR lines by Khomenko et al. (2003). This sup- 
ports the hypothesis of Sanchez Almeida & Lites (2000) and 
Socas-Navarro & Sanchez Almeida (2002, 2003): the study 
of the quiet Sun magnetism with only visible lines or only IR 
lines biases the resulting field strengths; the IN fields must be 
studied combining Stokes profiles from both spectral ranges. 

The inferred magnetic fields occupy 1.5% of the total 
FOV. Other observations with bett e r spat i al resolution (1"; 
'Socas-Nava rro & Sanchez Almeidal l2002t iKhomenko et alJ 
2003) yield a total filling factor of 1 %. We obtain a larger 
value thanks to the combined use of the information from vis- 
ible and IR lines, while all previous works used lines from a 
single spectral range. The filling fact or inferred for the same 
observations used in this paper by ISanchez Almeida et alJ 
(^003a) was only 0.85%. The difference is due to the dif- 
ferent model atmospheres used to reproduce the observa- 
tions. While here we u se a realistic multicomponent model, 
ISanchez Almeida et alJ (1200 3 a^ carried out Milne-Eddington 
inversions with a single magnetic componen t. Th us, the filling 
factors in pixels with mixed polarities (see § I5.2> are underes- 
timated. 

We have measured an average flux density over the FOV of 
9.6 G (those pixels without polarizatio n sig nals contribute to 
this average with zero field strength; § I5.1> . This flux is sim- 
ilar to that obtained by analysis of visible lines with 1" spa- 
tial resolution (see lSocas-Navarro & Sanchez Almeida 20021 
Lites 2002) while it is two times larger than the magnetic flux 
derived from IR lines ( Khomenko et al. 2003). Albeit our ob- 
servations have lower resolution (1 "5), the use of visible and 
IR lines reveals more magnetic structures than the use of lines 
from a single spectral range. The flux density is also larger 
than the magnetic flux obtained in § |3] using the magneto- 
graph equation applied to a single line. A combination of 
factors explain such bias. First, we obtained magnetic field 
strengths larger than the saturation limit of the lines, thus the 
circular polarization signals are no longer proportional to the 
magnetic flux and equation (Q is a lower limit of the real flux. 
Second, the presence of mixed polarities (see § I5.2> produces 
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cancellation of polarization signals and the magnetic flux is 
underestimated. Third, the lines weaken in strong magnetic 
concentrations, leading to a well known diminishing of polar- 
ization signals (e.g., Harvev & Livinaston 1969). Finally, the 
visible and IR signals often trace different concentrations co- 
existing in the resolution element. Then, combined inversions 
allow us to detect more flux than the inversion of one of the 
spectral ranges. 

The inversion indicates th e pr esence of mixed polarities in 
at least 20% of the FOV (§|52}- This is the largest fraction 
of mixed polarities obtained from direct measurements. The 
fraction of pixels with mixed magnetic polarities grows as the 
polarization signals weaken (Fig. lilt . More than half of the 
profiles with Stokes V signal lower than 0.2% of the contin- 
uum intensity have mixed polarities. 

We have estimated the probability of finding a given mag- 
netic field strength in the quiet Sun, i.e., the probabi lity d en- 
sity function of the magnetic field strength (PDF; § 15. 3> . It 
is largest for the weakest fields, however, an extended kG tail 
goes up to a limit around 1800 G. Despite the fact that the kG 
fields have only 30% of the total filling factor, they carry 75% 
of the total flux and 90% of the total magnetic energy. 

There is a clear imbalance of the magnetic flux towards 
one polarity. We measure a signed flux density of - 3.5 G 
which is mostly due to the strong kG fields (§ 15 .4> . The 
signed PDF ("Fig.ll3> shows how the weak fields are well bal- 
anced whereas the strong fields have a clear preference for the 
negative polarity. This imbalanced flux has the same sign as 
the sur rounding network, as it was already observed by Lites 
yOOZ). This imbalance is qualitatively similar to those im- 
posed in some MHD simulations of granular convection (see 
[Vogler 2003; Vogler & Schiissler 2003). These simulations 
have an initial unipolar field. The convection reprocesses this 
field to yield a continuous PDF with a clear imbalance in the 
strongest fields. 

The MISMA inversion of the IR profiles alone also gives 
rise to a PDF with extended kG tail (§ I6.4> . Such tail does not 
appear in Milne-Eddington inversions of the same profiles, 
which tend to select the component of weakest field in the 
resolution element. 

The magnetic fields studied in this paper fill only 1.5% 
of the total FOV. The unsigned flux density that we mea- 



sure is only a lower limit of the real flux. Larger fluxes 
are obtained in observations with bet t er spatial resolu- 
tion (iDominguez Cerdeiia et alJ l2003aHJ; ISanchez Almeidal 
2603b| and, furthermore, measurements based on the Hanle 
effect are compatible with a flux density o f 60 G o r more 
( Sanchez Almeida et aL. 2003b: Truiillo Bueno et al.l l200l 
Bommier et al. 2005). A large fraction of the magnetic 
structures remains undetected in Zeeman based observations. 
Based on our biased observations, we have produced a PDF 
with the distribution of magnetic field strengths. Obviously, 
this is not the real solar PDF. A natural step forward would be 
working out an unbiased PDF with the help of different mea- 
surements based on the Zeeman effect, the Hanle effect, and 
also MHD simulation s. This work of synthesis is carried out 
in a separate paper (Domingu ez Cerdeiia et a l. 2006). De- 
spite these and other efforts, we are still far from having a 
complete picture of the quiet Sun magnetism. There is a re- 
duced number of observations and each one is gathered with 
different instruments under different conditions. Thus, we 
need a survey of polarimetric data across many IN regions of 
the Sun, if we want to study the quiet Sun with a systematic 
and global point of view. Such data would allow us to start ex- 
ploring the connection of the IN with active regions, the chro- 
mosphere and corona and, eventually, to follow the evolution 
of the IN with the solar cycle. Observations with better spa- 
tial and temporal resolution are mandatory. They are coming 
up thanks to the advent of ne w generation gr ound based and 
space-borne solar telescope ([ Keil & Avakvaa.2003.) . These 
improvements, together with the development of more realis- 
tic numerical simulations, guarantee a rapid evolution of the 
field. 
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